Phase conversion in a weakly first-order quark-hadron transition 
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We investigate the process of phase conversion in a thermally-driven weakly first-order quark- 
hadron transition. This scenario is physically appealing even if the nature of this transition in 
equilibrium proves to be a smooth crossover for vanishing baryonic chemical potential. We construct 
an effective potential by combining the equation of state obtained within Lattice QCD for the 
partonic sector with that of a gas of resonances in the hadronic phase, and present numerical results 
on bubble profiles, nucleation rates and time evolution, including the effects from reheating on 
the dynamics for different expansion scenarios. Our findings confirm the standard picture of a 
cosmological first-order transition, in which the process of phase conversion is entirely dominated by 
nucleation, also in the case of a weakly first-order transition. On the other hand, we show that, even 
for expansion rates much lower than those expected in high-energy heavy ion collisions, nucleation 
is very unlikely, indicating that the main mechanism of phase conversion is spinodal decomposition. 
Our results are compared to those obtained for a strongly first-order transition, as the one provided 
by the MIT bag model. 



I. INTRODUCTION 



It is widely accepted that experiments in high-energy 
heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) have produced clear signals that nuclear 
matter undergoes a phase transition to a deconfined par- 
tonic phase at sufficiently large values of energy density 
[2, [2] . A similar transition presumably took place in the 
early universe a few seconds after the big bang [!, H| . In 
fact, this whole picture is expected from quantum chro- 
modynamics (QCD), which exhibits the phenomenon of 
asymptotic freedom. The nature of the quark-hadron 
transition, nevertheless, remains an open question. Al- 
though Lattice QCD Q, whose recently improved tech- 
niques allow for performing calculations with almost re- 
alistic quark masses @, seems to indicate a crossover, the 
possibility of a weakly first-order transition is not ruled 
out from the experimental point of view. In reality, most 
hydrodynamic calculations within high-energy heavy ion 
collisions adopt an equation of state which provides a 
strongly first-order transition [l| . 

Another point that is seldom mentioned is that, be- 
ing built on equilibrium assumptions, Lattice QCD ther- 
modynamics does not provide any information on the 
dynamical nature of the deconfining transition. In ac- 
tual experiments, the critical (dynamical) behavior could 
be very different from what one would expect from a 
crossover in the (equilibrium) phase diagram, simply be- 
cause the phase conversion is achieved via a nonequi- 
librium evolution process @, H[- Indeed, results from 
simulations in statistical mechanics show that the crit- 
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ical behavior can differ significantly if one compares the 
equilibrium phase diagram to the (noncquilibrium) time 
evolution of a given system, and suggest that the dy- 
namics could change the nature of the phase transition, 
even though the situation is still unclear 0, [lo| . For 
these reasons, the scenario of a weakly first-order decon- 
fining transition is physically appealing, especially since 
it comes out naturally by matching two equations of state 
that are realistic in their own regimes of temperature, 
namely the equations of state provided by Lattice QCD 
and by a hadron resonance gas. 

In this paper we investigate the process of phase con- 
version in a thermally-driven weakly first-order quark- 
hadron transition. We build an effective potential by 
combining the equation of state obtained from lattice 
simulations for one heavy and two light flavors of quarks, 
which we use for the partonic sector, with the equation of 
state of a gas of resonances for the hadronic phase. Us- 
ing standard techniques for the dynamics of first-order 
transitions, we compute bubble nucleation features, such 
as bubble profiles, critical radii, the surface tension and 
the free energy as functions of the temperature. We study 
the process of phase conversion evaluating the nucleation 
rate and investigating the time evolution of the temper- 
ature of the system and its hadronic fraction, as well as 
the role played by reheating. 

Bubble nucleation is one of different simplified mecha- 
nisms used to describe the dynamics of a first-order phase 
transition Q. In this kind of transition, for temper- 
atures slightly lower than the critical temperature, T c , 
the thermodynamic potential exhibits a metastable min- 
imum besides the global minimum. The former gradually 
disappears as the system cools down, and is turned into 
a point of inflection known as spinodal instability. By 
decreasing further the temperature, the transition fol- 
lows a qualitatively different (explosive) evolution, since 
the free energy barrier disappears and there is no extra 
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surface cost for small amplitude long-wavelength fluctu- 
ations. Therefore, depending on the dynamics of super- 
cooling, the phase conversion can proceed via more than 
one mechanism. In a slowly expanding system, the phase 
transition occurs through the nucleation of bubbles of 
the "true vacuum" state via thermal activation. In gen- 
eral, in order to know whether the system reaches the 
spinodal instability before nucleation is completed, it is 
necessary to investigate the time scale for thermal nu- 
cleation relative to that for the expansion. In the case 
of the quark-gluon plasma (QGP) presumably formed in 
ultrarelativistic heavy ion collisions, the expansion of the 
plasma is faster than that in the case of the early uni- 
verse during the cosmological quark-hadron transition by 
a factor of ~ 10 18 . Thus, it is clear that one should ex- 
pect major differences in the process of phase conversion 
when comparing the big bang to the little bang (llj . 

Our findings confirm the standard picture of a cos- 
mological first-order transition, in which the process of 
phase conversion is entirely dominated by nucleation, 
also for the case of a weakly first-order transition. On 
the other hand, we show that, even for expansion rates 
much lower than those expected in high-energy heavy ion 
collisions, nucleation is very unlikely, indicating that the 
main mechanism of phase conversion is spinodal decom- 
position. Our results are compared to those obtained for 
a strongly first-order transition, as the one provided by 
the MIT bag model. 

The dynamics of nucleation and spinodal decompo- 
sition in the hadronization of an expanding QGP af- 
ter a high-energy heavy ion collision has been studied 
under different approaches for more than twenty years 
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the case for the cosmological quark-hadron transition 
[23l |24| . However, previous studies were mostly based 
on equations of state obtained in the frame of the MIT 
bag model or within effective models, such as the linear 
sigma model [2^], the Nambu-Jona-Lasinio model [26| 
and the Polyakov loop model 27] in the environment of a 
strongly first-order transition. Only more recently, stim- 
ulated by the findings of Lattice QCD, a few studies have 
considered the case of a smooth crossover, and how its dy- 
namics compares to a strongly first-order transition, in- 
cluding the effects from fluctuations and inhomogeneities 
and the presence of a critical point [H, [29|, H(| . Our pro- 
posal in this paper is, as described above, to investigate 
a scenario which is in-between these two extrema. Fur- 
thermore, in the case of high-energy heavy ion collisions, 
in which there is a very fast quench down to the spinodal 
instability, the dynamics can possibly be well described 
by spinodal decomposition even in the case of a smooth 
(but very fast) crossover. 

The paper is organized as follows. Section II reviews 
the theoretical methods we use to describe the nucleation 
process of an expanding system, including a brief presen- 
tation of the thin- wall approximation following a different 
reasoning, and discusses the dynamics in an expanding 
background plus the role of reheating. Our numerical re- 



sults are discussed in Section III, where we display the 
relevant bubble features, as well as the time evolution of 
the temperature and of the hadronic fraction, and ana- 
lyze the process of reheating. Finally, our conclusions are 
presented in Section IV. 



II. THEORETICAL FRAMEWORK 

Thermodynamically, a phase transition happens when 
a given system shifts its state of equilibrium from one free 
energy minimum to another in response to the change of 
some critical thermodynamic parameter. In real transi- 
tions, this shift between equilibrium states is often an 
essentially non-equilibrium process and, in principle, one 
has no hope to describe it using the machinery of equi- 
librium thermodynamics. However, several natural sys- 
tems exhibit phase transitions in which the system is 
trapped in a metastable state (false vacuum) for a long 
time before reaching the final equilibrium configuration 
(true vacuum). This is certainly true in a strongly first- 
order scenario in which the time scale of supercooling 
(superheating) in a thermally-driven transition is not 
much faster than the typical reaction time of the sys- 
tem. A first-order phase transition occurs between non- 
contiguous states in the thermodynamic configuration 
space, and manifests itself as a discontinuity in the en- 
tropy This latent heat is a consequence of an energy 
barrier that prevents the system to simply roll down to 
the true minimum. If this barrier is high (characteristic of 
a strongly first-order transition), a statistically improb- 
able fluctuation is required, and the system is held in a 
metastable state for an appreciable time interval. 

A similar metastable behaviour sets in when the 
change in the critical parameter is fast as compared to 
the relaxation time of the system. This is the case of 
the hadronization of the QGP formed in high-energy 
heavy ion collisions, where the time scale of the expand- 
ing plasma is orders of magnitude shorter than the typi- 
cal reaction time of the deconfined matter. It turns out 
that it can be impossible (and immaterial) to determine 
the "real nature" of the phase transition in this dynam- 
ical system since the quenched plasma does not have to 
follow a continuous path in the thermodynamic configu- 
ration space. Assuming a weakly first-order scenario for 
the quark-hadron transition, we review in the next sec- 
tions some theoretical tools used in the description of the 
decay of metastable states. 



A. Homogeneous nucleation 

The paramount example of metastability in a first- 
order transition is the supersaturated vapor in a gas- 
liquid mixture below the condensation temperature. The 
mixture is composed of bubbles of liquid suspended in 
the gas, the two phases separated by a definite surface. 
Taking into account the surface energy of the partially 
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developed phase, the notion of thermodynamic equilib- 
rium can be extended in spite of the fact that the gas 
can not be in equilibrium with a fully developed liquid 
phase. Since it takes energy to make a surface, small 
bubbles will not be energetically favorable, and will be 
constantly created and evaporated. On the other hand, 
if a sufficiently large bubble is formed in the medium, it 
will tend to grow, eventually converting the entire system 
to the new phase. 

The mechanism of bubble nucleation can play an im- 
portant role in a first-order phase transition. It ignores 
the initial stages of the development of bubbles, being 
an effective theory for semimacroscopic elements of vol- 
ume. Still, one may have a flavor of what really happens 
by studying the dynamics of the nucleated bubbles. The 
question of the formation of bubbles is extremely rele- 
vant, mainly when external agents play the role of nu- 
cleating centers, leading to a significant increase of the 
nucleation rate. This is actually what happens to most 
natural systems, and is known as heterogeneous nucle- 
ation. However, we will not consider this kind of nu- 
cleation mechanism, but one in which bubbles originate 
from intrinsic thermodynamic fluctuations: the mecha- 
nism of homogeneous nucleation. It is a more fundamen- 
tal process for which there is a field theory which captures 
its basic features. Most of the formalism is based on a 
series of papers written by Langer in the late sixties [3l| , 
where the basic theoretical apparatus to describe the de- 
cay of a metastable state of a classical system interacting 
with a heat bath at temperature T is proposed (see also 
[32|). In this approach, stable and metastable phases ap- 
pear as local minima of a smooth energy functional E. 
Based on a phenomcnological droplet model, it is conjec- 
tured that, in going from one minimum to a neighbor one, 
the system is likely to pass across a saddle point which is 
a minimum of E in all directions of functional space but 
one, the latter giving rise to the instability. This saddle 
point configuration plays the role of the critical bubble in 
the formalism. Langer derived an equation of motion for 
the probability distribution of the system, and obtained 
a steady-state solution flowing across the saddle point, 
allowing for the calculation of a classical nucleation rate 
(probability per unit time per unit volume that a critical 
bubble nucleates) 
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where F is a prescribed analytic continuation of the free 
energy of the stable phase in the unbroken (single-well) 
version of the theory which becomes metastable after the 
analytic continuation. The parameter A_i is the negative 
eigenvalue of the fluctuation operator characterizing the 
instability of the saddle point configuration. 

The zero-temperature quantum field theoretic case, 
where transitions are exclusively due to quantum tun- 
neling, was considered by Coleman and Callan in (33|, 
and, when both quantum and thermal fluctuations act 
together, the decay rate is approximately given by (fT]), 



as proved by Affleck [34j | in a quantum mechanical con- 
text. Of course, the existence of a single direction of 
instability simplifies the extension of Langer's formalism 
to quantum field theories, and direct applications were 
proposed in the literature [HI, [36| . 

A good description of the nucleation process relies on a 
suitable choice of the free energy functional governing the 
dynamics of the bubbles. The natural choice corresponds 
to deriving the free energy from a more fundamental the- 
ory. When this is not feasible, and this is the case for the 
full QCD lagrangian, one can resort to a phenomenolog- 
ical approach, imposing symmetry requirements, as will 
be discussed in the next section. 



B. Effective potential and equations of state 

As customary, one can obtain information about the 
phase transition by studying the evolution of a scalar field 
<f> which represents the order parameter. It is reasonable 
to assume spherical symmetry for nucleating bubbles, so 
that one defines a coarse-grained free energy functional 
of the form: 



F{4>) = 47T J 



r 2 dr 



(2) 



Thus, the field </> evolves in space in the presence of an 
effective potential that can be parametrized in the form 
of a Landau expansion around the equilibrium phases, 



V(4>) = a{T)4> 2 - bTcj? + c4> 4 



(3) 



The parameters b > 0, c > and a(T), and the in- 
terpretation of the order parameter are determined by 
the scenario under consideration, namely the supercooled 
(T < T c ) QGP. One should notice that, in this simple ap- 
proach, the temperature enters only as a parameter of the 
effective potential. 

The potential is suitable for a first-order phase 
transition due to the properties of its extrema. The order 
parameter configurations (bubble profiles) are solutions 
of the following Euler-Lagrange equation: 
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(4) 



The potential ([3]) has two minima, <p q and 4>h, which 
correspond to the equilibrium phases. Here, 4> q (4>h) cor- 
responds to the quark (hadron) phase. There is a barrier 
separating <\> q and 4>h which can be associated with a la- 
tent heat, a jump in entropy from one phase to the other. 
In addition, one can prove that, as required by Langer's 
formalism, this theory has a saddle-point solution, 4>b, 
connectin g th e two minima and with a single unstable 
direction [3l|. The minimum <p q is conveniently chosen 
to be zero, while the other one is located at 



i (3bT + ^% 2 T 2 - 32a(T)c 



(5) 
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The transition temperature is reached when the pressure 
of the competing phases coincide. This condition sets a 
connection with thermodynamics through the identifica- 
tion 



p q (T)-p h (T) = V(<j> h ,T) ; V(0 q ,T)=O 



(6) 



where p q and ph are the pressures of the quark (decon- 
fined) phase and hadron phase, respectively. 

In order to proceed with the study of the phase conver- 
sion from the deconfined state to hadrons, one has to fix 
the potential, either by connecting it to pressures com- 
puted for each phase, as described above, or by extracting 
the effective potential directly from some effective field 
theory, such as the linear sigma model as, for instance, 
in Ref. [13, HH . We choose to follow the first procedure 
and, instead of using the MIT bag model as usually done 
for simplicity, we profit from the currently more robust 
knowledge of the equation of state of QCD in the two dif- 
ferent regimes (partonic and hadronic) which justifies the 
use of more realistic expressions for the pressure. Con- 
cretely, we use Lattice QCD results for Nf — 2 + 1 quark 
flavors to describe the high-temperature sector [33], and 
a gas of over 250 free resonances for the hadronic phase 
[38j |. This yields a weakly first-order deconfining tran- 
sition, to be contrasted to the usual case of a strongly 
first-order transition as provided by the bag model, and 
the value of the critical temperature is automatically de- 
termined by the crossing of the high and low tempera- 
ture pieces of the equation of state (EoS), as illustrated 
in Fig. 1. For comparison, we also use the bag model 
for the quark phase in our calculations and discussion, 
choosing the bag constant according to the critical tem- 
perature obtained by the crossing of the pressure curves 
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FIG. 1: Pressures for the hadron resonance gas (solid line), 
bag model (dashed-dotted line), and Nf — 2 + 1 Lattice QCD 
(dashed line) equations of state. 



Finally, the relations ([6|) leave two remaining condi- 
tions to be imposed on the potential. These can be cho- 
sen to fit a pair of physical properties of the critical bub- 



ble near the transition temperature, such as the surface 
tension and the correlation length [13| . 



C. The thin- wall approximation revisited 

It is always convenient to consider the following auxil- 
iary mechanical problem in the solution of Eq. ([4]): the 
one-dimensional dynamics of a particle moving in a dis- 
sipative medium in the presence of the potential (— V). 
Ignoring the dissipative term, strictly at T = T c , the so- 
lution of Eq. ^ can be written in terms of elementary 
functions. The aim of this section is to discuss approxi- 
mate analytic solutions for temperatures close (in some 
precise sense) to T c , where dissipation must be present. 

In order to explore this case, it is convenient to write 
the potential V as the sum of the critical potential plus 
a linear term which introduces deviations from T = T c 
[20l | . This is always possible for a quartic potential up to 
a shift in the zero of the energy. In the present case, we 
obtain V(4>) = W{4> + n) + Wq, where Wq is a constant, 



w{<j>) = c{4> 2 ~4>1) 2 +^ 



bT 

c 



a(T) 



a(T)bT 1 
~~c~ 2 4 



6T 

c 



(7) 

(8) 
(9) 



and n = —bT/Ac. The corresponding Euler-Lagrange 
equation has the same form as ([3]), with W' replacing 
V. The critical (T = T c ) potential (for which j = 0) is 
a symmetric double well with minima at cf> = +4>q. In 
terms of these parameters, the solution in the thin-wall 
approximation at T — T c is a kink which interpolates 
between the two symmetric minima: 



4> b (r) = o tanh 



R 



(10) 



where £ = l/((f)oV2c) can be thought of as a correlation 
length. When dissipation is neglected, the previous func- 
tion is a solution for any value of R. Given the specific 
form of the dissipation term, approximate solutions for 
the full potential are obtained for large values of 
The limit £/R <C 1 characterizes the so-called thin- wall 
regime [8j. 

For T < T c , the parameter j is negative and the two 
minima of W are shifted to 



h = V 



1 



J 



and 



(11) 
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When \j / (8c0q)| <C 1, the potential is only slightly tilted. 
In this range, the dissipation term will be small, and it is 
reasonable to assume that the interpolating solution still 
has a kink-like profile, i.e.: 



fair) = (j>f + 



1 — tanh 



R 



(13) 



where £ = yj2/c(4>t — 4>f) 1 - Still in this limit, one can 
show that the free energy of the bubble 4>b is given by 



F{<j> b ) 



AW + 4irR 2 a 



(14) 



where AW — W((/)f) — W(4> t ) « 2o|j'|, and the surface 
tension, a, is related to the parameters of the problem 
through 



dr 



dr 



3c£ 3 



(15) 



As discussed in the previous subsection, the value of a 
at T = T c and the constant value of £ will be considered 
as inputs in our treatment. In terms of these quantities, 
the coefficients of the potential ([3]) are completely deter- 
mined. Imposing that <f>b be stationary, we obtain the 
critical radius 



R,= 



2cr 

AW 



(16) 



Strictly speaking, the function (|13p is an exact solu- 
tion only at T = T c , where j = and R c diverges. A 
different approach consists in taking the dissipation into 
account in the following manner: (i) the solution con- 
serves energy; (ii) it is exact at T c . We can impose both 
conditions to a function of the form 



ipb(r) = 4>f + 7 
Condition (ii) gives 
Rc. = 



1 — tanh 



r - R c 



-i + 4c(M + 7 ) [0§-(M + 7) 2 ] 
whereas condition (i) is equivalent to: 



(17) 



(18) 



2 fdip b (r) 



dr = W{<p f ) -W(tp b (0)) . (19) 



r \ dr 

In the thin- wall limit (£/ ' R c <C 1), condition (i) is simply 
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w{4> f ) - w{4> f + 2 1 ) 



(20) 



Coupling (i) with (ii) leads to 7 = ((p v — cf)f)/2, £ = £ and 
R c = R c . Numerical results point to a kink- like solution 
with a decreasing escape value when one decreases the 



temperature. This fact explains the behavior of the sur- 
face tension plotted in Fig. 4 within this approximation. 

Once we have an (approximated) expression for the 
bubble solution, we can use it to calculate the decay rate 
shown in Eq. (fT]). This rate can be written as [3l| 



r = 



Vo _ 



(27T) 



AF/T 



(21) 



where AF is the difference in free energy between the 
saddle-point configuration (<fib) and the metastable phase 
(4>f). The prefactor To is a product of a dynamical fac- 
tor k (related to the expansion rate of the bubbles) and a 
statistical factor flo, which accounts for the first correc- 
tions to AF due to quadratic fluctuations around each 
extremum. As usual, these quadratic fluctuations are 
formally written in terms of the determinant of the fluc- 
tuation operator [V 2 — V" ((/>)], where <f> is either 0/ or 
<f>b- Thus, one obtains: 



n = (nRl t 



1/2 ^ 4irR 2 c a 



3/2 f det'[V 2 
\ det[V 2 - 



V"(<t> f )} J 

(22) 

where the prime indicates that the zero and negative 
modes are excluded. The first term in the r.h.s. of 
([22]) comes from the negative eigenvalue of the fluctu- 
ation operator along the direction of instability. The sec- 
ond term is the contribution from the three zero modes 
which are present because the bubble breaks transla- 
tional symmetry. Except for the four unpaired eigen- 
values of [V 2 — V" (</)/)], all the other delocalized eigen- 
values of that operator cancel the corresponding ones of 
[V 2 — V"(<frb))- Taking all the remaining contributions 
into account, it is possible to show that [35| 



1/2 



n = J-P) 3/2 



(23) 



where £ 9 is identified as the correlation length of the 
metastable phase. For the dynamical coefficient, Cser- 
nai and Kapusta derived the following expression (l3j : 



4g(C + 477/3) 
{AuoYRl 



(24) 



In the previous formula, £ and 77 are respectively the bulk 
and shear viscosity coefficients of the quark phase, and 
Aw is difference in the enthalpy density of quark and 
hadron phases (see also Ref. [4fj] for a careful discussion 
of the dynamical prefactor). 



D. Dynamics of the phase conversion in an 
expanding background 

As soon as a QGP is formed after a high-energy col- 
lision of heavy ions or in the early universe, it expands 
towards the empty space around it. Due to this expan- 
sion, the energy density and, therefore, the temperature 
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of the plasma drops and eventually becomes smaller than 
the critical temperature for the quark-hadron phase tran- 
sition. Assuming a weakly first-order transition, as dis- 
cussed before, the system then becomes metastable and 
the nucleation of bubbles of the cold phase is possible. If 
the expansion rate is large enough, however, the system 
supercools so fast that it may reach a thermodynamically 
unstable region of the phase diagram, the spinodal re- 
gion, before nucleation is able to drive most of the system 
to the true vacuum. In this case, the phase conversion 
is dominated by the process of spinodal decomposition, 
which is rather different from nucleation, since it is char- 
acterized by long-ranged fluctuations, instead of localized 
nucleating bubbles [8|. Here we focus on the nucleation 
stage during the expansion, and how the nucleation and 
expansion time scales compare. 

For simplicity, we consider a homogeneous and 
isotropic expansion of the fluid, and assume that the 
Hubble parameter H(x) = d^u^{x) is actually a con- 
stant. We also assume that the entropy is approximately 
conserved (globally) during the expansion, once the vis- 
cosity of the QGP is ver y lo w, as indicated by heavy-ion 
collision experiments (3, l4l|. However, the viscosity can 
not be strictly zero, once it is necessary for the process 
of nucleation at vanishing chemical potential (see, e.g., 
ref. [13|). For definiteness, we assume that the viscosity 
coefficient in eq. (|2"4"1) . £ + 4^/3, is of the order of the 
lower bound set by holographic models [42j ■ 

To have a measure of the importance of nucleation in 
the hadronization process, we estimate how much of an 
expanding plasma should hadronize through nucleation 
before the spinodal temperature is reached. In a first 
approach, we employ some simplifying assumptions. The 
first of them is to neglect reheating effects (which will 
be considered later) and, assuming the conservation of 
entropy: 



(25) 



The assumption of isotropic expansion leads to the Hub- 
ble law for the scale factor a(t), 



da , . y , 



(26) 



and Eq. (|26[) leads to an exponential decay of the tem- 
perature with time 



T{t) = T c exp(-Ht) , 



(27) 



where t = is chosen to correspond to T = T c . 

Given an equation of state and an effective potential, 
we can find the spinodal temperature of the system, T sp , 
and the time t sp the plasma takes to reach this temper- 
ature. For our equation of state, T sp /T c = 0.984, which 
leads to Ht sp ~ 0.01. The fraction f(t) of space which 



suffers nucleation from t(T c ) 
be overestimated as follows: 



until t sp = t(T sp ) can 



< 



4tt 



{t sp T c f ~ 10" 3 « 1 , 



(28) 



where R(t' , t) is the radius at time t of a bubble "born" 
at time t' [R(t',t) < c(t — t')]. In this estimate we also 
use the value H^ 1 = 10 fm/c, which should also give 
an extra overestimating contribution to f(t sp ) l . We see 
that, if we neglect reheating effects, typical QGP expan- 
sion times in heavy ion collisions rule out nucleation and 
therefore most of the plasma should hadronize via spin- 
odal decomposition, a result that was also obtained using 
the linear sigma model [lg, l2fj . 

However, once a bubble of the true vacuum is nucle- 
ated, an amount of latent heat proportional to the vol- 
ume of the bubble is released in the medium, so that the 
temperature does not fall exponentially as in the pre- 
vious case. If the nucleation rate is high enough, the 
released latent heat may win the competition against the 
energy loss due to the fluid expansion and the plasma 
reheats. This reheating can drive the system to temper- 
atures close to T c , decreasing the supercooling rate and 
considerably delaying its arrival at the spinodal temper- 
ature. In this case, when we may say that reheating is 
effective, the whole system is hadronized via nucleation 
of bubbles and T sp will be reached only some time after 
the transition is completed. 



E. Reheating 

In order to account for reheating effects, we make some 
assumptions aside from those cited on the previous sec- 
tions. First, the latent heat released in the formation of 
a true vacuum bubble is uniformly distributed through- 
out the whole plasma, which is consistent with that of 
bubble growth via weak deflagration [23] . Therefore, if 
s q (T) (sh(T)) is the entropy in the QGP (hadron gas) 
phase and s(T) is the space average of the entropy den- 
sity, 



s = s h f + (1 - f)s q , 
then entropy conservation implies 



o(0) 
a(t) 



s q (T c ) , 



(29) 



(30) 



where we set T — T c and s = s q at t — 0. For a 
supercooling 5 = (T c - T)/(T C - T sp ) < 1, we have 
s q (T) s q (T c ) and using the thermodynamic relation 
— p = AV — e — Ts, we find a relation between the tem- 
perature and the scale factor at a given time t: 



my 

T c J 



a(0)Y + Aa(T) 



*(*) 



s q (T c ) 



(31) 



f(ts P ) 



dt—R 3 (t,t sp )T[T(t)} 



1 H 1 is usually estimated to be in the range 1 — 10 fm/c for the 
expanding QGP created in heavy-ion collisions |4.'J . 
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where As(T) = s q (T) — Sh{T) is proportional to the la- 
tent heat I at T ps T c . Once we consider a constant 
Hubble parameter H, the scale factor grows exponen- 
tially, a(t) — a(Q)e Ht , and Eq. (j3Tj) becomes an implicit 
equation for the temperature as a function of time. We 
see that the first term on the r.h.s. of Eq. (|3Tj) accounts 
for the cooling due to the expansion, while the second is 
proportional to the hadronized fraction / and reflects the 
homogeneous reheating caused by the phase conversion. 
We compute / using the expression [4J| 



f(t) = 1-exp. 



' ; / ait 1 ) \ " „ r „. ... 47r „ . , 
fj£ I v / 1 i-ttv+'m i >■>!/■> 



(32) 

where R(t', t) is the radius of a bubble created at time t' 
with critical size R c [T(t')] at time t: 

R(t',t) = R c [T(t')]^ ) +£dt''v w [T(t'')]^- ) . (33) 

We suppose that the velocity of the bubble wall at the 
temperature T is given by the following Allen-Calm equa- 
tion, which relates the velocity of a domain wall to the 
local curvature H Wal : 



v w {T) 



Ap{T) _ AV{T) 



'1 



(34) 



very similar results for temperatures close to T c and be- 
come more and more different as the temperature is low- 
ered, as expected. One of the main differences is the 
broadening of the numerical profile, which does not hap- 
pen with the thin- wall result, indicating not only the fail- 
ure of this approximation for lower temperatures but also 
the failure of the nucleation picture itself as the system 
approaches the spinodal temperature. 



Numerical (T=0.995T c ) 
Thin-wall (T=0.995T c ) 
Numerical (T=0.9997T c ) 
Thin-wall (T=0.9997T ) 




5 10 15 20 

Radial coordinate [fm] 



25 



FIG. 2: Order parameter distributions (bubble profiles), 
where the broken (liadronic) phase is inside the bubble (r < 
R c ) and the unbroken (QGP) outside (r > R c ). 



where Ap(T) = -AV(T) = V q (T) - V h (T) > is the 
pressure difference between the two phases at a given 
temperature T. This expression corresponds to a steady 
growth of the bubble in which the pressure difference 
between both sides of the bubble wall is balanced by a 
damping force which is proportional to the velocity of the 
wall. The friction coefficient r\ is given by 



7/ = ffTaiT) 



where 



a(T) 



ar 



(35) 



(36) 



is the surface tension as a function of the temperature, 
and f] is a number of order one. 

We solve the set of equations (f3i~l) - ((3"4"|) numerically, 
and the results are presented in the following section. 



III. RESULTS AND DISCUSSION 

A. Bubble features 

In what follows, we calculate the main physical at- 
tributes of critical bubbles both numerically and using 
the thin-wall approximation. Fig. 2 shows bubble pro- 
files, $(r), for various temperatures. It is clear that the 
numerical results and the thin-wall approximation give 



From the bubble profiles we obtain two important 
quantities: the bubble critical radius and the surface ten- 
sion. The critical radius as a function of the temperature 
is shown in Fig. 3. Notice that the critical radius (de- 
fined here as the distance from the bubble center in which 
= 0o/2) diverges at the critical temperature, as seen 
in Section III CI and also diverges at the spinodal tem- 
perature. Notice, however, that for these high values of 
supercooling, the very concept of a bubble does not make 
sense any longer, once the width of the wall is larger than 
the bubble radius. 



30 
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FIG. 3: Radius of the critical bubble as a function of T/T c . 
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The surface tension, cr(T), exhibits an interesting be- 
havior, as displayed in Fig. 4. From its value at T = T c , 
a(T c ) — 2MeV/fm 2 , downwards it starts to grow up to 
a maximum value and then decreases from this point 
down. Roughly speaking, the value of the surface ten- 
sion is the slope p(r) = (d(f>/dr) 2 times the width of the 
wall £(T), i.e., the region in which p(r) ^ 0. As the 
temperature is lowered from T c , the wall is progressively 
broadened, but p{r) does not change much. As a con- 
sequence, the surface tension increases. It is only when 
p(r) starts to decrease faster than the increase of £(T) 
that <j(T) decreases too. This happens for lower temper- 
atures because the order parameter deep into the bub- 
ble, 4>{r = 0), becomes smaller, as can be seen in Fig. 2 
and, therefore, p(r) also becomes smaller. Notice that the 
thin- wall approximation is not very sensitive to this com- 
petition between the width of the wall and the disconti- 
nuity of the order parameter across it because its bubble 
profiles always connect the two minima, having 4>h(T) 
as the decisive parameter, which decreases more slowly 
than the exact (numeric) ef>(r — 0) for lower tempera- 
tures. Therefore, in the thin-wall approximation, <r(T) 
is a monotonically decreasing function of T. This non- 
trivial behavior of cr(T) suggests that the temperature 
in which cr(T) reaches its maximum can be interpreted 
as a (generous) limit to the applicability of the thin-wall 
approximation. Following the previous criterion, one can 
say that the thin-wall approximation fails for supercool- 
ing higher than 5 = 0.1. 



face tension makes this formula a good approximation 
to the exact value obtained using ([2]), once the functions 
R C (T), Ap(T) and a{T) are known. The different results 
to AF, calculated using a simple thin-wall (where a is 
a constant), using Eq. (|37p (which we call T-dependent 
thin-wall) , and from Eq. ([2]), are shown in Fig. 5. In our 
following computations, we use Eq. (|37p. 
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FIG. 5: Change in free energy, AF, calculated in three differ- 
ent ways (see text for details). The thin- wall approximation 
leads to an overestimate of the nucleation rate. 
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FIG. 4: Surface tension as a function of T/T c . 

Our aim in calculating all these quantities is to evalu- 
ate the change in free energy AF due to the presence of 
the bubble, which is an essential ingredient for the nucle- 
ation rate T (see Eq. (|2T|) '). AF can be calculated either 
directly from {5} or, in a computationally faster way, by 
using the formula 



4-ir 

AF(T) = -jRl(T)Ap(T) 



Ana(T)R 2 c (T) 



(37) 



This expression resembles the thin- wall expression for 
AF, but here the temperature dependence of the sur- 



B. Time evolution 

Now that we have all the necessary elements, we may 
investigate the time evolution of the system. We first an- 
alyze how the temperature evolves with time if we con- 
sider the reheating backreaction effects on the expanding 
system. As a first example, we consider three different 
expansion rates: H^ 1 = 100 fm/c (fast), H^ 1 = 600 
fm/c (critical) and H^ 1 = 10 7 fm/c (slow). From Fig. 6 
we can see that in the first case the temperature drops 
so fast that the system reaches the spinodal temperature 
without having time to grow enough bubbles to reheat. 
This is an example of a quench into the spinodal region. 
As the expansion rate H is lowered, the growing bub- 
bles eventually have time to effectively release enough 
latent heat. Then, the cooling process is reversed and 
the temperature is raised to a value close to (and lower 
than) T c . In this way the phase transition is completed 
through bubble nucleation. When the whole plasma is 
hadronized, there is no more release of latent heat, the ex- 
pansion once again dominates and the temperature falls 
abruptly. For the realistic equation of state we adopt, 
we find that the lowest (critical) expansion rate H at 
which there is reheating corresponds approximately to 
H^ 1 — 600 fm/c. Finally, we also consider a very low 
expansion and verify that the phase conversion proceeds 
almost entirely at T = T c , which is equivalent to hav- 
ing an almost in equilibrium phase conversion. Notice, 
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however, that, if we take a closer look at the tempera- 
ture as a function of time, even for the slowest expansion 
rate there is a slight supercooling in the beginning of the 
phase transition (as can be seen in the zoom for early 
times displayed in Fig. 7). In any case, the slower the 
expansion rate, the closer to equilibrium the transition 
evolves. Furthermore, the product Ht sp is the same for 
all expansion rates. 




0.98 



0.05 0.1 0.15 0.2 0.25 0.3 
Ht 



FIG. 6: Temperature, in units of T c , as a function of time, in 
units of 1/H, for different expansion rates. 



Dissipation effects can generally affect significantly the 
dynamics of phase conversion, even in the case of an ex- 
plosive spinodal decomposition [22| . In our framework, 
viscosity enters as a parameter in the dynamical prefactor 
of the nucleation rate (see Eq. Recalling that the 

viscosity of the QGP is presumably very small [3,[4l[ and 
that, on the other hand, holographic models set a lower 
bound (rjAds/s = l/4-7r) [42}, we compare the dynamics 
of the phase conversion for two values of viscosity: rjAds 
and Sr/AdS- Fig. 8 shows that by increasing the viscosity 
from rjAds to SrjAdS the supercooling decreases approxi- 
mately by a factor of two, although the time to complete 
the transition remains unaffected. 

The effects of different values of latent heat on the dy- 
namics of the phase transition was also studied. In Fig. 
9, we show the time evolution of the temperature consid- 
ering the two equations of state discussed in Section Hi Bl 
and an expansion rate H^ 1 = 600 fm/c. 

In Fig. 9 we can see that the EoS provided by the bag 
model, which leads to a latent heat larger than the one 
produced by the mixed EoS we adopt, implies a different 
transition dynamics. In fact, the larger latent heat of 
the bag EoS leads to a faster reheating of the plasma, 
since each nucleated bubble releases much more latent 
heat. Thus, a faster reheating means that the system 
will spend more time close to the critical temperature 
and the average nucleation rate will be much smaller for 
higher values of latent heat. Therefore, as can be seen in 
the figure, the nucleation process will take much longer 
to be completed. 
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FIG. 7: Temperature, in units of T c , as a function of time, in 
units of 1/H, for different expansion rates at early times. 
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FIG. 8: Temperature, in units of T c , as a function of time, in 
units of 1/H for different values of viscosity with H -1 — 600 
fm/c. 



We can now focus on the hadronic fraction of the 
plasma, i.e., the fraction which has been hadronized via 
bubble nucleation up to time t, f(t). As can be eas- 
ily seen from Eq. (|32[) , f(t) depends very strongly on 
5=1 — T/T c , another measure of the supercooling. In 
fact, for a small S, when the thin-wall approximation is 
valid, 



r(T) = ^exp[-AF/T]«^exp 



167TCT 3 



(38) 



which depends very strongly on 8. This implies that 
f(t) will also grow exponentially as the temperature de- 
creases, but before reheating effects are manifest. Indeed, 
after reheating the fluid temperature is near T c and it is 
unlikely that new bubbles nucleate. This means that 
after reheating f(t) grows due to the expansion of the 
existing bubbles, i.e. in a much milder fashion than dur- 
ing the supercooling phase, as can be seen in Fig. 10. 
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FIG. 9: Temperature, in units of T c , as a function of time, in 
units of 1/H, for the mixed and the bag model equations of 
state with H' 1 = 600 fm/c. 



FIG. 11: T(t) and f(t) for H' 1 = 600 fm/c. Full lines are the 
actual T(t) and f(t), and the dashed line is a zoomed T(t) for 
comparison with f(t). 



We also notice that, for a slow expansion, the first stage 
(nucleation of new bubbles) takes place very early, and 
almost all of the dynamics of phase conversion is due to 
the expansion of the bubbles. 




Ht 

FIG. 10: Hadronized fraction of space as a function of time 
for H~ 1 = 600 fm/c (critical expansion) and for _ff _1 = 10 4 
fm/c (slow expansion). 

By examining both T(t) and fit) on the same plot 
(Fig. 11), one can see quite clearly the relation between 
the reheating and the growth of fit) on one hand and, on 
the other hand, that after the reheating the temperature 
drops only after the phase transition is completed, i.e., 
after the system reaches / = 1. 

If the plasma is quenched directly into the spinodal 
region and its hadronization begins, there will also be 
a release of latent heat, just like in nucleation, and the 
system may also reheat. A curious possibility is that of a 
double process phase conversion, in which the reheating 
takes the plasma back to some metastable temperature, 
in which nucleation dominates. As a result, part of the 
plasma is converted through spinodal decomposition and 



part by (inhomogeneous) nucleation. It would also be 
interesting to explore experimental consequences of such 
a hybrid transition process. 

IV. CONCLUSIONS AND OUTLOOK 

Using two different equations of state, a realistic 
matching of hadron gas of resonances and Lattice QCD 
for Nf = 2 + 1 on one side, and the bag model (matched 
onto a hadron gas in the low temperature sector) on the 
other, we calculated both static and dynamical features 
of homogeneous bubble nucleation in a weakly first-order 
quark-hadron transition scenario, which is physically ap- 
pealing if one takes into account results from the lattice 
and from experiments, as well as the nonequilibrium na- 
ture of the phase conversion process, especially in the 
case of high-energy heavy ion collisions. 

After setting up an effective potential, we obtained nu- 
merically bubble profiles, critical radii, the surface ten- 
sion and the free energy as functions of the temperature. 
We also compared our numerical results to those derived 
in the thin-wall approximation. We showed how this ap- 
proximation (valid only near the critical temperature) 
can be used to overestimate the nucleation rate for lower 
temperatures. This indicates that an adequate approach 
to the dynamics of nucleation away from T c should take 
into account exact bubble profiles <fi(r), which in gen- 
eral have to be calculated numerically. With these static 
quantities we were able to feed a model for the dynamics 
of bubble nucleation in a homogeneously and isotropi- 
cally expanding plasma, which is often adopted in stud- 
ies of the early universe and heavy- ion collisions. Within 
this model, the only essential difference between these 
two physical settings is the value of the inverse Hubble 
constant H , which is about 10 19 fm/c in the early Uni- 
verse during the quark-hadron transition era, and about 
10 fm/c in the case of heavy-ion experiments. We com- 
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puted the temperature and the fraction of hadronizcd 
plasma as functions of time for different expansion rates, 
H, in order to make a quantitative estimate of the im- 
portance of nucleation in (homogeneous and isotropic) 
expanding systems. 

In the scenario of a weakly first-order transition inves- 
tigated in this paper, it is clear that nucleation remains 
as the dominant mechanism of phase conversion in the 
early universe, as expected from the standard cosmologi- 
cal picture. Nevertheless, previous estimates for relevant 
time scales using the bag model equation of state, which 
yields a much stronger first-order transition, may differ 
by a factor of two when compared to results from a more 
realistic equation of state. For high-energy heavy ion 
collisions, where the plasma expands very quickly, the 



main mechanism for phase conversion must be greatly 
dominated by spinodal decomposition, which possibly 
has some effect on particle correlations and fluctuation 
[IlliEliSlaEliSSEliiliilHSIEl- This issue will be 
addressed in a future publication. 
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